What if Dark matter is Bosonic and self-interacting 
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Abstract 

Recently the problem of singular galactic cores and over-abundant formation 
of dwarf galaxies, inherent to the standard cold dark matter model, had at- 
tracted a great deal of attention. One scenario which may be free of these 
problems invokes a self-interacting Bose-field. We find the limiting core den- 
sity in this model due to the self-annihilation of the scalar field into its own 
relativistic quanta. The limiting density may correspond to the observable 
one if there is only one dark matter component. Alternatively, there may be 
more than one dark matter species and the annihilation of one species may 
be very efficient with subsequent expansion of the other, thus avoiding the 
problem of singular cores. 
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Introduction. The nature of the dark component which constitutes most of the matter 
in the Universe remains unknown. It is usually assumed that the dark matter (DM) is cold 
and non-interacting. Adding a non-zero cosmological constant to the non-baryonic dark 
matter gives rise to the so-called ACDM model. In spite of its great success in reproducing 
the observed Universe at large scales, in the recent few years it has become more and more 
evident that the ACDM model suffers from predicting too much power of density pertur- 
bations on small scales. First, the central density distributions in dwarf and low surface 
brightness galaxies seem to have finite cores in contrast [0 with the singular profiles pre- 
dicted by N-body simulations (see however Ref. ||). Secondly, high resolution numerical 
experiments predict that galaxies as our own should contain several thousands clumps of 
dark matter of the size of a dwarf galaxy Q . However, only a few of these are observed. 

Several solutions to this problem have been suggested in the literature: i) the dark matter 
is not cold, but warm [||; ii) the dark matter is interactive ||; Hi) the dark matter is warm 
and interactive J7J; iv) the power spectrum of ACDM has sharp drop on subgalactic scales 
||; v) the dark matter is in the form of a self-interacting scalar field 0. All these ideas are 
not new and have been previously considered in different applications. For instance, warm 
dark matter was considered in flO|| , self-interacting dark matter was studied in Ref. |TT 
while the scalar field model of Ref. PI was considered in Ref. fT2 



Clearly all these models require an appropriate degree of fine-tuning in their parameters 
(e.g. mass and self-couplings). For example, the dark matter component of Ref. || should 
interact sufficiently strongly to allow the dark matter in the halo to start relaxing on time- 
scale of the age of the Universe, but the relaxation should not be too fast and already 



completed otherwise a gravothermal instability will inevitably occur [13], leading to the 
formation of singular profiles. This is not a problem within the scalar field model where the 
pressure which supports the core is generated by the self-repulsion of dark matter particles. 
However, for the halo dominated by the scalar field the radius of its core is only a function 
of the field mass and self-coupling [|T^,[| and therefore they should obey certain relation to 
fit observations. 

In this paper we will study further constraints on the parameters of scalar field model 
with mass m and self-coupling A (with interaction potential defined as A0 4 /4) and give a 
testable prediction for the core radius and core density of the gravitational bound object 
composed by scalar quanta. In particular, we analyze constraints which are due to inevitable 
self-annihilation of the self-interacting scalar field without conserved quantum numbers into 
its own relativistic quanta. 

Self-interacting scalar field and gravitationally bound clumps. It is assumed that the 
scalar field is in the state of coherent oscillations, just like the familiar axion in cosmological 
setting, and that the gravitationally bound object is formed out of the field clump. Of 
the prime interest for comparison with observations is the case when self-gravity of the 
field dominates (generalization for the case of the scalar field in the external gravitational 
well, created e.g. by baryonic component, is straightforward). The parameters of such 
configuration can be estimated analyzing a simple equation of "hydrostatic equilibrium" in 
non-relativistic limit [[EJ 

dP(r) _ p(r)M(r) 
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= Wp(r) . (2) 

dr 

The pressure P(r) and the density p(r) have to be understood here as quantities averaged 
over the period of field oscillation, which gives for the equation of state |P~2"|,|9[ P ~ Ap 2 /3m 4 



These equations have a simple solution p(r) = p(0)sin(x)/x, where x = r/r c and the core 
radius satisfies^ 



m 4 M|j 



67rr^ 



(3) 



or 



r c ~ 10 2 >/A(eV/m) 2 kpc. (4) 

Note that the density in the core, p(0), is unrelated to the core radius, but the density is 
limited. Equation of state changes to relativistic, P = p/3 when O > m/vA, where O is 
the amplitude of field oscillations, and the object became unstable [[L2[]. This corresponds 
to the maximum density, p = m 2 0Q, 

Pmax ~ rn 4 / A . (5) 

This density corresponds to the object which is close to the limit of black hole formation. 
Since this is the only restriction on the density if the field has conserved quantum numbers 
(which corresponds to the complex scalar field), the model with conserved particle number 
is unappealing (the excess dark matter cusps may Bose-condense into black holes in the 
centers of galaxies fL4| , but this can not be a solution for all systems). Therefore, here we 
consider the model where the dark matter particle number is not conserved. For example, 
the real scalar field can decay or self-annihilate. The corresponding constraints on the 
density arising from self-annihilation, which are inevitable because of the non-zero A, will 
be considered below. Let us first though consider under which conditions the coherent field 
configuration forms. 

Condensation. If the self-coupling is not very large, during the collapse of spatially 
large configurations the coherence of the field will be destroyed even if the configuration 
develops from the small overdensity in initially homogeneous oscillating field. Namely, after 
virialization at each spatial point there will be many streams of particles, each with different 
vector of velocity. In such situation the analysis above and Eq. (||) are not applicable. If 
velocities dominate, the virialized configuration will behave in the gravitational field as non- 



interacting usual dark matter. However, with time the Bose-condensation will occur |L4 
This will happen even with very small values of the self-coupling since the scattering is Bose- 
enhanced if the phase-space density of particles is larger than unity, which will be always the 
case with parameters satisfying Eq. (|]). If the relaxation time for the Bose-condensation is 
smaller than the age of the Universe, we may apply Eq. (01) for the final configuration. 



1 This definition of r c coincides with Ref. 12] but is y/E smaller compared to Ref. ||, because in 
P] the r c was defined via p(r c ) = 0. 
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Let us consider now the relaxation in virialized clumps due to the scattering process 

-l 



20 following Refs. [14],[15]. The inverse relaxation time is t 



R 



[1 + n)apv e m , where 

a is the corresponding cross section, and v e is the characteristic velocity which characterizes 
the depth of the gravitational well. We are taking into account the possibility that mean 
phase-space density of 4> particles, n, can be large, which accounts for the factor (1 + n) 
in the expression for the relaxation time. For particles bound in a gravitational well, it is 
convenient to rewrite this expression in the form H 



t 



R 



rrr 



1 + 



m A vl 



(6) 



Let us consider first the case n > 1. The relaxation time in the clump is smaller than a 
given time r if 



P > 



m 



(7) 



With p = O.O2M /pc 3 ~ 5 x eV 4 and v e = 100 km s _1 , which approximately corre- 
sponds to the parameters of the cores of dwarf galaxies, the relaxation time will be smaller 
than the age of the Universe if 



A > 10 



-15 



m/eV) 



7/2 



(8) 



If this condition is satisfied, the Bose condensate in the center of gravitational well should 
form [ I4"| . 



In addition to the condensation due to self-interaction, there is a process of purely grav- 
itational relaxation, also with subsequent formation of the coherent field configuration, see 
Refs. 
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However, this should be efficient only if density inhomogeneities are large (e.g. 
during initial stage of the collapse), while we are interested in continuing condensation in 
the already virialized clump as well. In addition, in all range of parameters relevant for the 
present discussion, the condensation solely due to self-coupling is efficient. 

Self-annihilation of the condensate. Besides taking into account the specific Bose- 
enhancement during the process of collisional relaxation, one has also to consider the decay 
of the condensate. This is also peculiar and may correspond to a "laser" effect flI2UI7| . In 
other words, the Bose-stimulation of relevant processes has to be taken into account. 

In such a case it is important to know up to which distances the decay or annihilation 
products stay in the resonance with each other. If the initial configuration is not condensed, 
but can be described by some distribution of (free) particles over momenta, one can employ 
the Boltzmann kinetic approach [I7||. In a stationary state (created, say, by the process of 
violent relaxation) the distribution of particles in the phase space is a function of integrals 
of motion like energy. This does not mean though that one can neglect the red-shift in 
the collision integral for particles which are at the same energy level. This is because 
particles are moving. However, there is enough phase-space available for the decay products 
to stay in the resonance if the distribution of "parent" particles is isotropic [17]]. With the 
average momentum of the distribution of the parent particles going to zero (condensation), 
the collision integral start to diverge. The Boltzmann approach breaks down and one can 
employ |1| the formalism of particle creation by time dependent classical background. 
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For the condensed state in equilibrium in the gravitational well we have the following pe- 
culiarity: the red-shift cancels out and the decay products stay in the resonance throughout 
the whole configuration. First, in the condensed state particles are not moving. Second, in 
the state of hydrostatic equilibrium the gravitational energy and internal interaction energy 
are tuned precisely in the way to cancel the gravitational redshift. This can be understood in 
several ways. Indeed, in hydrostatic equilibrium, the total mass of configuration M (which 
includes all forms of energy) takes its minimum possible value under perturbations of internal 
structure. Some small amount of particles can be moved around resulting in 5M = 0. As a 
Gedanken experiment one can think to move around particles which are going to decay. The 
condition 5M = means that at infinity the products of decay will have the same energy. 
Therefore, they stay in resonance on WKB trajectories labeled by quantum numbers at 
infinity. In other words, from the point of view of local observers, in hydrostatic equilibrium 
/iy^oo = const, where \i is chemical potential, see e.g. Jji|. Compare this with the energy of 
free moving particle as measured by the local observer, uo^fgw = const. Therefore, /i/u — 
const and the decay process stays in the resonance. 

Let us make this explicit for the self-annihilation of the Bose-condensed field in the metric 
ds 2 = fdt 2 — f~ 1 5ijdx l dx\ where / = 1 + 2$ and $ is Newtonian gravitational potential. 
Condensed state can be described by the field configuration of the form = <f)(t)ifj(x), where 
<p(t) is periodic function of time which to the first approximation is <p(t) = (j) smujt, with 
to close to m. Let us describe quanta which propagate in this metric as WKB plane waves, 
5<p oc g v {t)e %VjX \ where pj = Pj(x) are spatial momenta which satisfy g^ v p^p v = m 2 s or 
Po = / 2 |p| 2 + f m cs = const. Here m e ff includes contribution to the particle mass from the 
interaction with the condensate, mg ff = m 2 + 3A0o-?/> 2 /2. Mode functions satisfy the equation 

8'<j) - f 2 V 2 5(j) + f(m 2 + 3A0)50 = . (9) 

Therefore the equation for g p (t) is 

9 P + P 2 o9 P - 2cu 2 (l + 2$)ip 2 qcos(2ujt)g p = , (10) 

where q = 3A0q/4c<j 2 is the resonance parameter. We see that change in the gravitational 
parameter $ with the distance gives insignificant change in the "effective" resonance param- 
eter. More important is the change in ip 2 (x), but most of the modes with the same value of 
Po can stay in the resonance up to the distances comparable to the core size. The important 
fact which leads to this conclusion is that u = const. 

The problem reduces to particle creation by time dependent homogeneous classical back- 
ground and is very well studied, see e.g. JL9|. We briefly review the results below. The num- 



ber density of the created particles grows exponentially with time, = exp where the 
characteristic exponent is positive and non-zero in narrow resonance bands and its numerical 
value is model dependent, being a function of the coupling constants of the theory. There 
can be several channels of decay, e.g. can be coupled to photons leading to peculiar cosmic 
maser effect discussed in Jl2| which by itself puts constraints on the strength of the coupling 



to the electromagnetic field. The strength of the coupling to the electromagnetic field 
(as well as to any other possible but hypothetical field) is unknown and might be even zero. 
However, we cannot disregard the decays of the condensate into its own quanta because of 
the self-coupling A. 
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Let us consider the consequences of the condensate decays. The rate of particle pro- 
duction as a function of particle momenta k is determined by the growth rate of unstable 
solutions of the Mathieu equation for the corresponding mode functions 

g k + [A-2qcos(2T)}g k = 0, (11) 

and at the center of the iV-th instability band the parameter /i^r is given by 

N 

^ N = 2N{2 N - i {N-iy.y ' ( } 

The coupling of <fi to the electromagnetic field gives Jl2| A = 4k 2 /m 2 and q = 2k(p /mf l j ) . 



The decay of (j) into two photons is saturated in the first instability band which is centered at 
A — 1. In such a case, \i\ = qm/2 and the products of the decay have momentum k = m/2, 
with a width 5k « /i. 

For the self-annihilation 40 — > 2<j), as we've seen one finds A = (k 2 + m 2 )/m 2 + 2q and 

q = 3\<f) 2 /4m 2 . (13) 

The self-decay of occurs in the second instability band centered at A = 4 with /z 2 = q 2 m/16, 
momentum k = \^3m and width 5k ~ [i. Clearly, the products of the self-annihilation are 
ultra-relativistic and easily escape the gravitational well. 

The rate of decay p is a function of the amplitude of the field oscillations 0o and therefore 
is a function of the energy density in the core: p = m 2 ^. Notice that the exponential growth 
of the particle number in the resonance bands is due to the Bose-statistics: already created 
particles stimulate production of new quanta. However, particles which leave gravitational 
well do not participate in the stimulation process. Therefore, the exponential growth occur 
only if D = fir c > 1, where r c is the core radius of the field configuration. If D ^> 1 initially, 
the number density of decay quanta will grow exponentially in time in the region of the 
core [^] because in each decay process two identical particles are produced and travel in 
opposite directions. In a sense, this system is equivalent to the inversely populated laser 
medium placed between reflecting mirrors. The resulting explosion will reduce the density 
in the core below the level at which D ~ 1. If the core was growing gradually, starting from 
the small density, which is likely the case in astrophysical situation, the density will just stop 
growing at the condition D « 1 even if the infall of particles continues due to condensation 
from surrounding non-relativistic "gas" . In this regime the luminosity in relativistic particles 
will be equal to the rate of Bose-condensation. 

Let find the maximum core density corresponding to the condition D ~ 1 in the case of 
self-annihilation. The condition q 2 mr c ~ 1 with r c defined in Eq. (|J) and q defined in Eq. 
(0) givesQ 



2 Numerically, the assumption of perfect hydrostatic equilibrium turns out to be not vitally im- 
portant here. Had we used instead of r c the distance over which the gravitational redshift equals 
to the width of the resonance band, we would have obtained exactly the same relation, but with 
power 1/2 replaced by the power 2/5 in the prefactor of m /A. 



6 



(771 \ ^ 777^ ( 1 \ 777^ 
s^xj -~feJ -• (14) 

Interestingly, there are indications that the halo central density is nearly independent of 
the mass from the galactic to the galaxy cluster scales, with average value of around p = 



O.O2M /pc 3 0. With this value of density, Eq. (0) gives 



A « 1(T 8 (m/eV) 7/2 . (15) 

Comparing this with the condition (||), we see that the Bose-condensation is efficient indeed 
in this parameter range. 

Let us finally estimate the maximum core density which corresponds to non-stimulated 
self-annihilation (D < 1, and the rate of Bose-condensation is small). The core density will 
stop changing effectively when the rate of annihilation will became comparable to the age 
of the Universe, to- We find 

I 771 

P ^W^~- (16) 

We see that not only the core radius, but also the core density may be uniquely defined 
in terms of the mass and the coupling of the scalar field. This provides a severe test of the 
model. The core radius or the maximum density (|T4|) will be changed somewhat, if the field 
is embedded in an external gravitational well (created by baryonic matter), but this will not 
alter the required parameter range significantly. 

However, the picture may be more complicated and perhaps more interesting with core 
radius (^) and limiting density ( |14D not related to the observed characteristics of the dark 
halos. Assume that initially the field contributes to the most of the dark component 
and that the parameters are such that the Bose-condensation is efficient on a time scale 
shorter than the age of the Universe. Finally, suppose that in the core of the halo which 
forms the condition D ^> 1 is satisfied. This is a modification of the scenario considered 
in Refs. [0,|T7|,|T4[ where the possibility of electromagnetic radiation was advocated. Here 



annihilation of the field into itself will do the job, but we require that the Bose-condensation 
with subsequent inflow into the core and decays are efficient enough to re-process the major 
fraction of the dark halo into relativistic particles. The baryonic core (and any other dark 
component) will then expand after loosing part of gravitating central mass. The problem of 
singular cores may be avoided in this way in a wider range of parameters. 

Needless to say, our proposal of annihilating dark matter, as a solution of CDM problems 
on small scales, is not limited to a self-annihilating scalar field, but applies as well to other 
possible forms of dark matter particles which may annihilate into different species. 

We are grateful to V. A. Berezin for useful discussions. 
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